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We model repetitive quantum error correction (QEC) with the single-error-correcting five-qubit 
code on a network of individually-controlled qubits with always-on Ising couplings, using our pre¬ 
viously designed universal set of quantum gates based on sequences of shaped decoupling pulses. 
In addition to serving as accurate quantum gates, the sequences also provide dynamical decoupling 
(DD) of low-frequency phase noise. The simulation involves integrating unitary dynamics of six 
qubits over the duration of tens of thousands of control pulses, using classical stochastic phase noise 
as a source of decoherence. The combined DD/QEC protocol dramatically improves the coherence, 
with the QEC alone responsible for more than an order of magnitude infidelity reduction. 


I. INTRODUCTION 

Quantum error correction [1^4] (QEC) makes it possi¬ 
ble to perform large scale quantum computations with a 
finite error rate per qubit [SHTTj. Much like their classical 
counterparts, quantum error correcting codes (QECCs) 
rely on adding redundant qubits to control errors. How¬ 
ever, unlike, e.g., in the classical information transmis¬ 
sion problem, qubits remain subject to errors all the time, 
in particular, during the syndrome extraction. Thus, to 
achieve scalability, special fault-tolerant (FT) protocols 
must be used both for QEC and for any operation with 
the encoded qubits. This increases the overhead and is 
one of the reasons why the error probability thresholds to 
scalable quantum computation are so small—e.g., around 
1% per local gate in the case of the toric and related 
surface codes [5J HU [T3]. The number of qubits needed, 
measurement complexity, and stringent requirements for 
gate speed and fidelity are also among the reasons why 
an experimental demonstration of the repetitive quan¬ 
tum error correction with a universal quantum code so 
far remains an elusive goal. |14It25] 

A possible way to loosen these requirements is to com¬ 
bine active QEC with one of the passive error reduction 
techniques depending on the correlations in the dominant 
decoherence channel P5l434) . In particular, dynamical de¬ 
coupling (DD) is most effective [55MT] in suppressing the 
effects of low-frequency (e.g., 1//) noise which is often 
the leading mechanism for the loss of phase coherence. 
Moreover, DD can be used to achieve scalability in gate 
design, since pulses and sequences intended for a large 
system can be constructed to a given order in the Mag¬ 
nus series [42] on small qubit clusters |37l |38] . DD is also 
excellent in producing accurate control for systems where 
not all interactions are known as one can decouple in¬ 
teractions with the given svmmetrv |431 144] . and it can 
be implemented to remain stable when environment has 
sharp spectral features |15| or high-frequency modes |4S]. 
or even with substantial pulse errors |i71 H5] . In short, at 
least in principle, using DD at the lowest level of coher¬ 
ence protection could substantially reduce the required 
repetition rate of the QEC cycle. In many instances, this 


could make a crucial difference, enabling the use of QEC. 

Recently we made a substantial progress toward de¬ 
veloping a combined DD/QEC coherence protection pro¬ 
tocol by constructing a universal set of quantum gates 
based on soft-pulse DD sequences. iiaisD| The gates are 
designed to work on a network of qubits with always- 
on Ising couplings forming a sparse bipartite graph. In 
addition to providing accurate control, these gates also 
work as decoupling sequences, suppressing the effect of 
low-frequency phase noise to second order in the Mag¬ 
nus series. With these gates, we demonstrated [15] the ef¬ 
fectiveness of repetitive QEC using single-error-detecting 
QECC [[4,2,2]] encoding two qubits in four by simulat¬ 
ing full unitary dynamics of five driven qubits in an Ising 
chain, using low-frequency classical noise as the source of 
decoherence. 

We have also studied[5D] the errors associated with the 
gates similar to those constructed in Ref. HSj In a sys¬ 
tem with always-on pairwise qubit couplings, for any gate 
constructed perturbatively to a given order K, only the 
errors forming clusters that involve up to AT -I- 1 qubits 
are suppressed. Large-weight clusters of correlated errors 
can be suppressed exponentially when gates are executed 
fast enough. However, such a choice can only be made 
with a sufficiently sparse coupling network. Increasing 
the maximum degree z of the connectivity graph with 
pulse duration and other gates’ parameters fixed may 
lead to proliferation of large uncorrectable error clusters. 

In our previous workjl^, we simulated a linear Ising 
chain with z = 2, an arrangement most favorable for con¬ 
trolling multi-qubit correlated errors |5T]. On the other 
hand, the optimal arrangement for surface codes is pla¬ 
nar. The corresponding analytical bound on maximum 
gate duration needed for FT is small [50]. Thus, it re¬ 
mains an open question whether perturbation-theory- 
based gates like those constructed in Ref, can be prac¬ 
tical for use in repetitive QEC. 

In this work we simulate numerically repetitive quan¬ 
tum error correction using our universal gate set |491 [50] 
in a network with z > 4 as would be needed for the 
surface code. Specifically, we use a six-qubit star graph 
(see Fig. [^ of Ising couplings between the qubits with 
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FIG. 1. (Color online) Arrangement of qubits on a bipartite 
star graph for implementing the [[5,1,3]] code. This particular 
arrangement of the qubits was chosen to minimize the total 
number of gate operation during the stabilizer generator mea¬ 
surement cycle, (a) The qubit to be encoded, Qe, is initially 
placed at the center. At the end of encoding, it is swapped 
with an ancilla qubit, initially Qi. (b) During the stabilizer 
measurement cycle, the single ancilla qubit Qe at the center 
is used to sequentially measure all four stabilizer generators 
supported by the five qubits around it. 


z = 5, and simulate QEC with the [[5,1, 3]] code both in 
the traditional and in the error-detection (Zeno) modes. 
This code can actually be seen as a variant of a surface 
code with rotated periodicity vectors. [52l |53] We simu¬ 
late full unitary dynamics over several error correcting 
cycles (up to seventy thousand shaped pulses) with in¬ 
stantaneous ancilla-based projective measurements, and 
classical Gaussian correlated phase noise used as a source 
of decoherence. We consider the cases of low-frequency 
noise with Gaussian correlations, as well as a bimodal 
noise generated as a combination of a low-frequency and 
a high-frequency components. The constructed protocols 
show substantial improvement of coherence compared to 
the case of unprotected qubits, including over an order 
of magnitude average infidelity reduction attributable to 
error correction alone. 

The structure of the paper is as follows. In Sec. [TTj 
we give a brief overview of dynamical decoupling, of the 
universal gate set used, and our implementation of the 
[[5,1,3]] code on a star graph. In Sec. Ill we present the 
results of numerical simulations. We give Gonclusions in 
Sec.HVl 


II. GATE AND CODE IMPLEMENTATION 
A. Dynamical control on an Ising network 

The goal of dynamical control is to drive the desired 
unitary evolution of a quantum system over a given time 
interval. While the details of the dynamics during the 
interval may differ greatly, the net result of such an evo¬ 
lution can be, to some extent, insensitive to the details of 
system’s interactions. For example, in the simplest case 
of single-qubit dynamical decoupling, the qubit interac¬ 
tions are averaged out during the period of the controlled 
evolution. 


We consider the following general Hamiltonian 

H = Hc+Ho, Ho = Hs + Hb + Hsb. (1) 

where He = Hc{t) is the time-dependent control Hamil¬ 
tonian, and the remaining Hamiltonian Hq is separated 
into the parts Hg and Hb acting on the qubit “system” 
and on the bath respectively, and the system-bath cou¬ 
pling Hamiltonian H^b- 

In this work, following Refs. HHI and (SUl we consider 
a sparse bipartite network of qubits with the Ising cou¬ 
plings between nearest neighbors [M]. 

( 2 ) 


and decoherence due to slow dephasing of individual 
qubits, generally described by the following bath and 
bath-coupling Hamiltonians: 

HB='Y^Bi, Hsb = 

i i 

Each qubit is assumed to have its own individual bath, 
i.e., the bath operators Bj commute with each other, and 
the coupling operators Ai commute with each other and 
all Bj , j yf L 

The decoupling technique assumes that the control 
Hamiltonian Hq dominates the dynamics. We implicitly 
assume that any large energies have already been elimi¬ 
nated from the system Hg and system-bath coupling Hsb 
Hamiltonians by a rotating wave approximation (RWA). 
Then, the Hamiltonian ([^ can be viewed as an effec¬ 
tive Hamiltonian for any set of interactions as long as 
the transition frequencies of the neighboring qubits dif¬ 
fer sufficiently. Similarly, the bath model <§ is an effec¬ 
tive description of qubits operating well above the bath 
frequency cut-off to eliminate direct spin flip transitions, 
with dephasing caused, say, by phonon scattering. 

We also assume the ability to control the qubits indi¬ 
vidually, with the control Hamiltonian 

Hc^Y.Hc\ = l E (4) 

j fj-=x,y,z 

where the time dependent control signals Vj^{t) are ar¬ 
bitrary, except for some implicit limits on their ampli¬ 
tude and spectrum. Our gates [121110] are designed as se¬ 
quences of one-dimensional pulses, with the control fields 
on a given qubit applied along x, y, or z axis exclusively, 
so that only one of 14 (t), Vy{t), or 14 (t) can be non¬ 
zero at any given time. We also imposed a restriction 
that no pulses be applied simultaneously on any pair of 
neighboring qubits. 

As a result of these assumptions, the multi-qubit 
unitary evolution operator with the complete Hamilto¬ 
nian 0 over the duration of a single-pulse interval can 
be written as a product of mutually commuting terms. 
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each of them involving the controlled qubit and vari¬ 
ous products of Pauli operators for its uncontrolled 
neighbors. m Each of these operators can be computed 
order-by-order in the time-dependent perturbation the¬ 
ory; in Ref. [SD] we carried such an expansion up to 
third order. In each order of the series, the depen¬ 
dence on the pulse shape is encoded in terms of just 
a few coefficients [35H1Z1 fSD] . For example, the first- 
order correction is expressed in terms of just two such 
coefficients, the time averages of cos and sin 
over the duration of the pulse, where (j){t) = dt'V{t') 
is the time-dependent rotation angle corresponding to 
the given pulse shape V{t), 0 < t < Tp, and Tp is 
the pulse duration. With generic pulse shapes, such 
as a Gaussian, this produces errors that scale linearly 
with Tp. Specially designed self-refocusing pulses can 
be constructed to suppress this effect, e.g., to linear 
or quadratic orders in powers of Tp, depending on the 
shape [571 Wj\ . For example, to the linear order, this is 
done by choosing a functional form V(t) which guaran¬ 
tees {cos (j){t)) = {sin = 0. If the pulse shape is 
symmetric, V{t) = V{Tp — t), this requires only one ad¬ 
ditional condition on the shape [5711571155] . 

While in a multi-qubit setting such special pulse shapes 
do not eliminate all first- or second-order errors over 
the pulse duration, the resulting series have fewer terms 
which can be subsequently dealt with easier by properly 
designing the pulse sequences. 


B. Universal gate set 

With generic set of inter-qubit couplings, increasing 
the number of qubits requires progressively longer se¬ 
quences to decouple the inter-qubit couplings [35]. How¬ 
ever, when the couplings form a bipartite graph, such 
a decoupling to an arbitrary (hxed) order can be done 
with a single sequence of a finite duration independent of 
the number of qubits in the system[35]. We constructed 
a gate set formed by an arbitrary single-qubit rotation 
and an entangling controlled-phase gate [more precisely, 
arbitrary-angle ZZ rotation, exp(—iacrfcrl) for a pair of 
neighboring qubits i and j]. According to general the¬ 
ory, such a gate set is universal [57]. These gates can 
be executed simultaneously on an arbitrary set of non¬ 
neighboring qubits (pairs of qubits), and in addition pro¬ 
vide DD protection for every qubit. In particular, all the 
Hamiltonian terms not directly involved in the gate are 
averaged out. 

Single-qubit gates SiliO] are based on the leading- 
order dynamically corrected gates[5HJ [55], in turn based 
on the Eulerian path construction[60j . The original 
single-qubit DCG sequence[5Hl [55] guarantees leading- 
order cancellation of an arbitrary bath coupling. This 
is achieved by executing a sequence of identically-shaped 
pulses driving the qubit through an (extended) Eule¬ 
rian cycle on the Cayley graph corresponding to the 
decoupling group. Explicitly, the single-qubit DCG 


sequence!; 


can be formally written as 


(X)(/)(Y)(/)(X)(/)(Y) (Y)(X)(Y)(X) (P), (5) 

where (X) and (Y) represent finite-duration tt pulses in 
X and Y direction, (P) the pulse nominally implement¬ 
ing the desired single-qubit rotation, and / is the com¬ 
posite pulse implementing a unity operation by running 
a half-time double-amplitude version of P followed by an 
identical pulse applied in the opposite direction. As writ¬ 
ten, the sequence works for pulses of arbitrary symmetric 
shapes, as long as these shapes remain the same during 
the sequence. 

To build dynamically-protected single-qubit gates on 
a bipartite qubit network with always-on couplings, we 
separated the DCG sequence into a part to be executed 
on the sublattice A [X pulses in the original sequence 
([^] and a part to be executed on the sublattice B {Y 
pulses in the original sequence replaced by X pulses). 
Each of these are partial-group sequences as they go over 
Eulerian cycles corresponding to subgroups of the two- 
sublattice decoupling group, specifically chosen to control 
Ising bath coupling As a result, the entire sequence 
is only effective against de phas ing, and it requires self¬ 
refocusing pulses (see Sec. HA) to achieve the leading- 
order cancellation. 

The construction allows simultaneous rotations in any 
set of non-neighboring qubits (e.g., the entire sublattice 
A or B can be rotated at once), with P representing the 
desired rotation or zero applied field on idle qubits. In 
actual implementation we used the stretched pulse P of 
duration 2Tp, so that the unity operation I in Eq. (j^ is 
composed of two pulses of duration Tp. Overall, the du¬ 
ration of such a single-qubit gate is Ibr^. The Hadamard 
gate P[ is implemented as a product of two rotations, 
with the net duration 32rp. 

Same sequences used with second-order self-refocusing 
pulses [571137] (see the portion t < 16rp on Fig. yields 
second-order cancellation of inter-qubit couplings and the 
bath terms in Eq. (|^, except for terms proportional to 
the commutators [Bi,Ai\. These terms are readily inter¬ 
preted as the derivatives of time-dependent fields acting 
on the qubits. Such terms can also be canceled [50]. e.g., 
using symmetrized versions of our sequences (involving 
32 pulses instead of 16), if one uses more complicated 
pulse shapes akin to those developed in Ref.jOTJ While se¬ 
quences achieving higher cancellation orders can be read¬ 
ily designed using the same general approach [5H11501105] . 
the advantage of the particular sequences we use in this 
work is that they are shorter. 

Two-qubit ZZ-rotation gates[301]50] are designed 
using a different approach, see Fig.[^ The idea is to selec¬ 
tively decouple some of the inter-qubit interactions, with 
the needed rotations generated by the residual interac¬ 
tions when the sequence is repeated over some specified 
amount of time. This only requires conventional decou¬ 
pling sequences which are, generally, easier to design. 

The qubits are divided into four sets: idle qubits on 
sublattices A and B (depending on the chosen graph). 
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FIG. 2. (Color online) Pulse sequences used to implement the CNOT gate between qubits Qs and Qq on a star graph. It is 
a combination of four single-qubit gates (0 < f < 16rp, 32Tp < t < 48rp, 48Tp < t < 64Tp, and 64rp < t < 80Tp) and a single 
instance of the ZZ-coupling sequence, 16Tp < t < 32Tp. For better accuracy this latter sequence has to be repeated several 
times, we used N^ep = 5 repetitions, see Ref. [50] for detailed discussion of the associated errors. Second-order self-refocusing 
pulse shapes Qi(7r) and Qi{n/2) from Refs. (37] and |47| are used. The shading shows the direction of the applied pulses as 
indicated. 


and the sets A! Q A and B' B which together make 
up all of the pairs of neighboring qubits where we want to 
preserve the couplings. The corresponding sequences are 
denoted Va, Vb, Va', and Vb>- The universal idle-qubits 
sequences Va and Vb must decouple both the system Q 
and the bath ([^ Hamiltonians, and have sufficient flex¬ 
ibility so that the coupling with a neighboring opposite- 
sublattice qubit driven by the sequence Vb' and Va' re¬ 
spectively could also be decoupled. On the other hand, 
the sequences Va' and Vb' executed on the pairs of qubits 
to remain coupled must average out the bath Hamiltoni¬ 
ans ^ , but leave some fraction / of the original coupling 
([^ between these qubits. 

We designed the global sequences Va and Vb to allow 
for construction of local versions of the sequences Vb' (/) 
and VA'if), with some range of allowed fractions /. This 
makes the fraction / locally adjust able [SHI, to accommo¬ 
date for possible local variations of the couplings Jij. In 
this work we assume all couplings equal (non-zero Jij = J 
iff the sites i and j are connected), and use the fastest 
version of these sequences of duration Tgeq = 16rp with a 
fixed fraction / = 1/2, as used originally in Ref. [321 Over 
the duration of the sequence, for each pair of qubits des¬ 
ignated to be coupled, the original coupling Jij = J in 
Eq. (U) is reduced to /J = J/2, which gives the rotation 
angle a = /JTsec/2. 


We constructed a CNOT gate using the identity [551IM] 

= y.WjVjVi, exp(-^^a/a/)n (6) 

= ZjXkYkexp(^-i^a^a{jYk, (7) 

where, e.g., Xu = and Xk = respectively are 

the unitaries corresponding to ±7r rotations of the A:-th 
qubit around the X axis. Eq. Q requires a ZZ rotation 
with the rotation angle a = Tr/d. Thus, the pulse dura¬ 
tion Tp and the number of repetition A/ep must satisfy 
the design equation [35] 

A^rep JTp = — . (8) 

Larger values of A^rep correspond to smaller values of the 

perturbation-theory parameter Jt^ which improves the 
fidelity as it provides better decoupling. On the other 
hand, this also increases the cost in terms of the number 
of pulses. The actual set of driving helds used to imple¬ 
ment the CNOT gate with = 1 are shown in Figj^ 
For our calculations we used A^rep = 5. 

We also implemented two other controlled two-qubit 
gates using the identities 

= e-^^/^X^Z.Zu exp(-*%;a^)Xfc, (9) 

^ exp(-7^a|a|), (10) 
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as well as the swap gate as a sequence of three cnot 
gates. m 


C. Five-qubit code on a star graph 

We use the smallest single-error-correcting code 
formally denoted as [[5,1,3]]. This distance-three 
code encodes a single qubit in a two-dimensional sub¬ 
space Q of the 2®-dimensional Hilbert space of n = 5 
qubits. It is a stabilizer code [2]: the subspace 

Q = m : G,1V') = |V'), ^ = l,2,...r} 

is a common -|-1 eigenspace of the r = 4 independent 
commuting stabilizer generators. 


Gi = XZZXI = a^a^alal, 
G 2 = IXZZX = 

G 3 = XIX ZZ = alalalal, 
Gi = ZXIXZ = 


( 11 ) 


expressed as Kronecker products of single-qubit Pauli op¬ 
erators erf, /i = x,y,z. Notice that to reduce the confu¬ 
sion with the pulse unitaries in Sec. |II B here we quote 
both the commonly used positional and the traditional 
notations for multi-qubit Pauli operators. 

As for any stabilizer code, encoding of the five-qubit 
code can be done efficiently [2]. We have used the con¬ 
ceptual encoding circuit in Fig. [^a), which produces 
the code in the basis with the logical operators X = 
—XXXXX and Z = ZZZZZ. This circuit is based 
on a representation of the five-qubit code as a code word 
stabilized (CWS) code[ 68 ], and was constructed as a sim¬ 
plification of the circuit containing the Hadamard gate 
on the information qubit, encoder for the classical five- 
qubit repetition code, and the graph state encoder[Hi^. 
Explicitly, the resulting basis wavefunctions correspond¬ 
ing to the eigenvalues = (— 1 )™, m = 0 , 1 , are (up to 
a normalization) 


= jOOOOm) - jOllOm) -k JlOOlm) - Jllllm) 

-k jOOlOm) -k jOlOOm) - JllOlm) - llOllm) 

X (lOOOlm) -H 11110m) -H jOlllm) -H JlOOOm) 

-k 10011m) - jOlOlm) - JlOlOm) -f JllOOm)). (12) 

To implement the same circuit on the star graph, we 
used two more SWAP operations, plus an additional SWAP 
at the end to place the ancilla at the center, see Fig.j^b). 
This initializes for the stabilizer generator measurement 
cycle shown in Fig. [^ 



| 0 > - 

| 0 ) -{h} 

| 0 ) -{h} 

| 0 ) -{h} 

| 0 ) ^ 

(b) lb) - 

FIG. 3. (a) Conceptual encoding circuit for the [[5,1,3]] code 
using the Hadamard, CNOT, and controlled-phase gates. On 
the input, the first four qubits are initialized in jO) states 
and the last qubit contains the state \ip) = ajO) -I- djl) to 
be encoded. On the output, the corresponding logical state 
a 10) -I-/31 i) of the five-qubit code is produced. The decoding is 
done by inverting the encoding circuit, (b) Actual encoding 
circuit implemented on a six-qubit star graph. Two SWAP 
operations are required since qubits on the leaves can only 
interact with the qubit 6 in the center. In addition, the ancilla 
qubit is SWAPped to the center at the end to prepare for the 
measurement cycle, see Fig. [^ 
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pulse sequences described in Sec. El and classical zero- 
mean Gaussian phase noise with Gaussian correlations, 

(Mt)) = 0 , {Ait)A,it')) = (13) 

as a source of decoherence, cf. Eq. Notice that for 
a single uncontrolled qubit, such a field would produce 
asymptotic dephasing rate I/T 2 = (•00/2) cr^T„. 

The corresponding many-body unitary dynamics has 
been simulated with a C++ program using the EigenS 
library [TH] for matrix algebra. The program uses a 
custom-built algorithm to schedule the pulse sequences 
and measurement events, and the fourth-order Runge- 
Kutta algorithm for explicitly integrating the time de¬ 
pendent Schrodinger equation for the unitary time evo¬ 
lution of clusters of multiple qubits. In all simulations 
shown we used 1024 time steps per nominal pulse du¬ 
ration Tp, resulting in relative integration errors better 
than 10“®, comparable to numerical precision. 


III. SIMULATIONS 

We implemented the described encoding/decoding and 
the measurement circuits at the Hamiltonian level, using 


A. Quantum error detection mode 

We first consider the working of the [[5,1,3]] code in 
the error detection mode (quantum Zeno cycle [TTl 17^ 1. 
In an actual experiment, one is supposed to measure the 
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FIG. 4. The circuit for a single cycle of measuring the stabilizer generators (111 for the hve-qubit code on a star graph. On 
the input and the output, the first hve qubits contain an encoded state. The same ancilla qubit 6 at the center is used for each 
measurement; M stands for measurement and resetting to |0) if needed. We implemented this circuit which uses 16 CNOT gates 
(duration 144rp each) and 8 rounds of Hadamard gates (duration 32Tp each) applied in parallel, with the total measurement 
cycle duration of 2560rp (640rp per stabilizer generator). In our simulations, the entire cycle is repeated several times for 
repetitive QEC. 


stabilizer generators repeatedly, with the experiment ter¬ 
minated once an error is detected as indicated by a non¬ 
zero syndrome bit. In our simulations, instead, each syn¬ 
drome measurement was replaced by an instantaneous 
projection 


Po = (1 + ^6')/2, (14) 

which selects the many-body sector with the ancilla qubit 
Qq at the center in the state |0). The success probability 
averaged over the initial state \tj}) was calculated accord¬ 
ing to the expression 

Po = Tr{UpoU^Po), (15) 


where U is the N x N unitary evolution matrix up to the 
moment of measurement, po = M~^Pm is the density 
matrix describing the uniform distribution of the initial 
wavefunctions in a subspace of dimensionality M, and 
Pm is the corresponding projector. In our simulations, 
iV = 64 is the dimensionality of the six-qubit Hilbert 
space, and we compute a reduced NxM evolution matrix 
V which include only M = 2 columns corresponding to 
the number of basis states of the initial qubit, see the 
encoding circuit in Fig.|^ Respectively, we used Eq. (15) 
in the form 


Po = ^Tr{V^PoV). 


(16) 


Given the reduced evolution matrix V = V{t) &t a given 
time moment t, and the corresponding ideal evolution 
matrix Vo, the overall fidelity averaged over the initial 
state can be calculated using the expression 

f = + (17) 

“ ^ M{M + \) ■ ^ ’ 

The derivation is similar to that given in the Appendix 
of Ref. 27] for the case of M = N. 


Results of simulations for several sets of parameters 
of the Gaussian noise, the r.m.s. amplitude a and the 
correlation time t„, are shown in Fig.[^ Each plot is an 
average over 20 instances of the stochastic noise, with the 
time axis starting at the first measurement after the end 
of the encoding. Having in mind an experiment where 
the success probability and the state fidelity would be 
measured separately, and to match the quantities com¬ 
puted in Ref. |49l we plot the success probability (SP), 
and the fidelity “with measurements” (WM) conditioned 
on the error-free syndrome measurements, Psucc = F/po, 
where po is given by Eq. (16). Notice that this expression 
is an approximation which ignores possible correlations 
between po and Psucc- These correlations would be ab¬ 
sent with ideal syndrome measurements; we expect them 
to be small in our case since the measurement fidelity 
is high. The effect of such correlations is additionally 
suppressed since Psucc is numerically close to one. 

To compare the contributions of the DD protection and 
of the projective measurements (Zeno cycle) to the over¬ 
all fidelity, in Fig. we also show the average hdelity 0 
calculated when decoupling pulses are applied but “no 
measurements” are done (NM), and when no projective 
measurements and “no pulses” are applied (NP). Since 
they involve no projective measurements, these quanti¬ 
ties are independent of the success probability (15). For 
each version of the cycle, filled symbols show the infideli¬ 
ties 1 — F after each syndrome measurement, while open 
symbols show the corresponding infidelities to the end of 
the final decoding. Notice that thus computed fidelities 
involve all six qubits; final infidelities could be addition¬ 
ally reduced by tracing out all but one information qubit, 
see Sec. IHIBI 

These plots show about an order of magnitude infi¬ 
delity reduction due to QEC during the cycle. The code 
can detect any one- and two-qubit error, and a small 
fraction of higher-weight errors. The fact that the Zeno 
cycle works, indicates that errors seen by the code are 
not dominated by multi-qubit correlations. In addition. 
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FIG. 5. (color online) Infidelities during the Zeno cycle for 
different noise correlation times and noise amplitudes as in¬ 
dicated. The time axis starts after the completion of the 
encoding circuit, at the instance of the first measurement, see 
Figs. Sb) and 1^ The different curves correspond to cases 
where no pulses are applied (NP), DD pulses are applied but 
no measurements are done (NM), and with the syndrome mea¬ 
surements (WM). Closed and open symbols respectively rep¬ 
resent the infidelities at the end of each syndrome measure¬ 
ment and at the end of the final decoding. Note that the axis 
for the cumulative success probability (SP) is on the right. 


the infidelities increase sharply with shorter noise corre¬ 
lation times, as expected due to the asymmetry of the 
single-qubit gates, see Sec. II B[ 

Two of the plots shown have exactly the same noise 
parameters and use the same pulse shapes as in our ear¬ 
lier work[^ where Zeno cycle was simulated with the 
[[4, 2,2]] error-detecting code, with five qubits arranged 
in a chain. The corresponding success probabilities and 
state fidelities are similar in magnitude. We believe this 
to be a combined result of an improvement due to more 
efficient code and faster syndrome measurements in the 
present case, negated by increased errors due to larger 
connectivity of the star graph, as discussed in detail in 

Ref. 1501 


B. QEC mode 

In this mode we simulated projective measurements of 
the ancilla by applying instantaneous projection opera¬ 
tors Po [Eq. ( [I4| )] or Pi = 1 —Pq. These are six-qubit pro¬ 
jectors selecting the sector with the ancilla qubit Qe at 
the center in the state |0) or |1), respectively. Given the 
normalized wavefunction 41 of the system, the projectors 


should be chosen with the probabilities po = ('I'| Pq |4>) 
and Pi = (4'|Pi |4>) = 1 — po, respectively. This im¬ 
plies a separate simulation would be needed for every 
state of the initial qubit (see the encoding circuit in 
Fig.§. Instead, to speed up the simulations, we calcu¬ 
lated the reduced unitary evolution matrix V and used 
the probability (161 averaged over the initial state of the 
qubit. The normalization of V was corrected after each 
projection. This approximation is similar to that used 
in the previous section to define the fidelity Psucc condi¬ 
tioned on the string of zero-syndrome measurements in 
each previous cycle. Here, we also expect the effect of 
any potential unaccounted correlations to be suppressed 
due to the smallness of pi = 1 — po- 


As in the previous section, we simulated decoherence 
with classical phase noise applied on all six qubits in¬ 
volved in the simulations. The noise was uncorrelated 
between different qubits. For each qubit, the noise was 
generated as a stationary zero-mean Gaussian stochas¬ 
tic process with Gaussian time correlations. Individual 
traces of such simulations for three realizations of the 
Gaussian stochastic noise with identical correlation time 
T„ = 32Tp and different r.m.s. amplitudes as indicated 
are shown in Fig. Each panel shows four different 
infidelity measures 1 — P computed during a single sim¬ 
ulation run. The fidelities Fh and are computed ac¬ 
cording to Eq. 0, respectively, just before and right 
after each projective measurement. The fidelities and 
computed at the same time moments include idealized 
recovery channel, see Eq. (181 and the discussion below. 


The five-qubit code is a “perfect” single-error- 
correcting code, since the fifteen (2'^ — 1) non-zero syn¬ 
dromes corresponding to four stabilizer generators © 
are in a one-to-one correspondence with the fifteen single¬ 
qubit errors. We used this idealized map for decoding. 
Notice, however, that in our simulations the stabilizer 
generators are measured sequentially, with the entire 
measurement cycle typically repeated just a few times. 
To increase the syndrome measurement fidelity, we did 
not adhere to a fixed measurement cycle and instead trig¬ 
gered the beginning of a cycle by a syndrome measure¬ 
ment returning a non-zero bit. After the fourth mea¬ 
surement, the correction would be computed and applied 
immediately. Typically, the infidelities 1 — Pb computed 
right before a trigger event were small, whereas right af¬ 
ter the infidelity 1 — Pa jumps to near one, as the wave- 
function is projected outside of the code. The infidelities 
1—Pb remain large right before the subsequent three mea¬ 
surements, creating an easy to spot four-dot pedestal in 
the combined trace. For exam]de, no trigger events hap¬ 
pened in the top trace in Fig. ^ (cr = 5 x 10“^/rp), one 
happened in the middle trace [a = 20 x 10“^/rp), and 
two in the bottom trace (cr = 50 x 10“^/Tp). 

To look beyond the simple system fidelity ( [T7| ), we also 
calculated the fidelity including the idealized recovery 
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FIG. 6. Sample error correction traces for [[5,1,3]] code in the 
presence of the stochastic phase noise on all six qubits. The 
noise correlation time is r„ = 32rp and noise amplitudes o 
are as indicated (in units Tp^). Plots show different infidelity 
measures 1 — F during the stabilizer measurement cycle and 
at the end of the decoding. Here F\, and Fk are the regular 
fidelities [Eq. ( |17[ )], respectively, computed just before and 
right after each projective measurement. Fj( and F^ are the 
corresponding fidelities which include idealized recovery for 
single qubit errors, see Eq. (181. Larger symbols indicate the 


infidelities at the end of the decoding c ircu it, where the circle 
corresponds to the full-system fidelity (171, and the diamond 


to the single-qubit fidelity after tracing out the qubits away 
from the center. 


map. 


15 


F'{V, Ho) = F{V, Vo) + Y, E^Vo), (18) 


where FCj, i = 1,...,15, are all single-qubit errors on 
the peripheral qubits, and F{V,Vo) is the usual fidelity 
(17). This expression corresponds to idealized error cor¬ 


rection, with the summation over all single-qubit errors 
corresponding to that over all possible syndromes. 

We should mention that in our calculations both fi¬ 
delity expressions include the ancilla qubit which has not 
been traced out. However, since the ancilla is reset to |0) 
state after each projective measurement, it is effectively 
excluded for the fidelities Fa and F^ computed right after 
the measurement. The ancilla is also included in the full- 
system fidelity F computed at the end of the decoding 
circuit, but not in the final fidelity F' which only looks 
at the state of the single qubit in the center. In our plots 
these fidelities are shown with bigger symbols. 

The plots in Fig. [^show the simulated infidelities for 
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FIG. 7. Infidelities 1—F for the [[5,1, 3]] code averaged over 25 
realizations of the stochastic Gaussian noise, with the noise 
correlation times r„ and amplitudes a as indicated. Large 
infidelities after trigger events have been excluded from the 
averages. F are the regular fidelities [Eq. 0] computed right 
after each projective measurement, F' are the corresponding 
fidelities which include idealized recovery for single qubit er¬ 


rors, see Eq. 


In addition Fn and F'd, respectively, are 


the fidelities ( |17P and (18) for DD-only simulations with the 
same pulse sequences run but error correction turned off (no 
projective measurements). Larger symbols indicate the re¬ 
spective quantities at the end of the decoding circuits, with 
F' and F{, replaced by the average single-qubit decoded fi¬ 
delities, with the qubits away from the center traced out. 


one simulation run each, they are strongly affected by the 
details of the particular noise realization and the mea¬ 
surement results simulated probabilistically. In Fig. we 
show (in the logarithmic scale) infidelities averaged over 
25 different realizations of the stochastic noise. To re¬ 
duce the unphysical fluctuations, large infidelities after 
the trigger events have been excluded from the averages. 


The data in Fig. also include average infidelities 
1 — Fd and 1 —F)^ [Eqs. ( [T7| ) and (18)] produced in identi¬ 
cal simulations but with error correction turned off (the 
same pulse sequences but no projective measurement). 
Except for the plots in Fig.j^d), where QEC becomes rel¬ 
atively ineffective due to strong noise with shorter corre¬ 
lation time, the DD-only infidelities show a substantially 
steeper growth than those where both DD and QEC was 
active. The overall QEC effectiveness can be quantified 
by comparing the final single-qubit infidelities 1 — F'j^ 
and 1 — F' at the end of the decoding (two larger cir¬ 
cles). The corresponding ratios of average infidelities for 
different panels in Fig. [^are: (a) 14.3, (b) 15.0, (c) 9.73, 
and (d) 1.36. Except for the data in Fig. Hd), QEC in 
these plots gives average infidelity reduction by an order 
of magnitude or better. Notice that for this data, trigger 
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events are rare; here QEC fidelity is similar to that for 


the Zeno cycle, see Sec. Ill A 


In the three plots where QEC works well, Fig.[^a)-(c), 
the data for \ — F' is some two orders of magnitude below 
that for 1 —E, indicating that in the present setup single¬ 
qubit errors strongly dominate. This is in an apparent 
contrast with the results of our Ref. [SOI where we con¬ 
cluded that multi-qubit errors are an unavoidable feature 
of the perturbatively designed gates. We notice, how¬ 
ever, that due to asymmetry of single-qubit gates, in the 
presence of time-dependent noise, the leading-order error 
terms are single-qubit Pauli operators [SD], with the coef¬ 
ficients scaling as a derivative of the classical fields Ai(t). 
Further correlated errors are formed in higher orders of 
the Magnus series, they can be represented as connected 
clusters on the connectivity graph. On the star graph, 
these include a single bond joining the ancilla at the cen¬ 
ter to one of the code qubits, and, in the next order, two 
bonds, which could result in a correlated error involving 
the ancilla and two qubits of the code. Thus, after the 
ancilla is projected during the measurement, the remain¬ 
ing errors on the qubits forming the code are expected to 
have smaller weight than they would with a different con¬ 
nectivity graph. The applicability of these arguments is 
improved by our choice Ai^ep = 5, which gives the pertur¬ 
bation theory parameter Jr^ec = where Tgec = 
is the typical sequence duration, see Eq. Q. 


This analysis is confirmed by the plots in Fig.|^ which 
show infidelity traces for different amplitudes of the noise 
with the correlation time = I28rp. The two top pan¬ 
els, with noise amplitudes cr = 0 and tr = 5 x I0“^/rp, 
show near identical plots for F', indicating that with the 
noise parameters as in Fig. [^b), multi-qubit errors are 
strongly dominated by the systematic errors due to the 
couplings J. At the same time, single-qubit errors are 
dominated by the stochastic noise, since at tr = 0, the 
plots for F and F' fall nearly on top of each other. 


Similar conclusions can be also drawn from the data 
in Fig. 1^ which shows the effect of a much faster noise, 
with the correlation time r„ = r^. Namely, the infideli¬ 
ties in Fig. I^a) were generated by averaging the results 
of 25 simulations with different realizations of Gaussian 
noise with the correlation time = 128rp, while the 
noise for infidelities in Fig. ib) , in addition, also in¬ 
cluded weaker but faster-varying noise components with 
Tn = Tp. Dynamical decoupling has nearly no effect on 
such a fast noise. Respectively, the usual infidelity 1 — F 
increased by an order of magnitude, while the inhdelity 
I — F' including idealized recovery map [see Eq. (18)] in¬ 
creased by more than two orders of magnitude. Such a 
different scaling of the two infidelities dominated by one- 
and multi-qubit errors, respectively, is consistent with 
the expectation of the absence of DD protection against 
the faster noise. Quantitatively, the ratios of the final 
average single-qubit infidelities at the end of the decod¬ 
ing in runs with and without error correction are 15.7 in 
Fig.|^a) and 8.0 in Fig.j^Jb). 



X h = 128 , o = 20 X 10"^ 
n p 
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FIG. 8. As in Fig. but with longer noise correlation time, 
Tn = 128rp. Different infidelities are labeled as in Fig. 


T /t = 128 , a = 20x 10" 
n p 



x^/Tp = [1,128],a = [2,20]x10"' 



FIG. 9. As in Fig. comparing the effect of fast dephasing 
noise. The plots at the left show infidelities averaged over 
25 realizations of a Gaussian stochastic noise with r.m.s. am¬ 
plitude (T = 20 X 10“^/rp and r„ = 128rp. For the data 
on the right, in addition, there was also a weak noise with 
CT = 2 X 10“^/Tp and a much shorter correlation time = Tp. 
Such a noise is not affected by the DD. With the addition of 
the fast noise, the infidelity 1 — F increased by about an order 
of magnitude, while the infidelity 1 — F' accounting by multi¬ 
qubit errors increased by two orders of magnitude, consistent 
with the expected absence of DD protection against the fast 
noise component. 
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IV. DISCUSSION 

For many years, the road to building a quantum com¬ 
puter appeared to be straightforward: one just had to 
manufacture a sufficient number of quality qubits and 
implement a universal set of quantum gates of sufficiently 
high fidelity. Now that we are there, or nearly there, it 
turns out that fidelity is not the ultimate measure of per¬ 
formance in large qubit systems. As the number of qubits 
in a quantum computer grows, exponentially so does the 
number of ways it can go wrong. To understand what 
is going on in a particular implementation of a quantum 
computer would require detailed numerical simulations, 
including as many qubits and as much physical detail as 
possible. 

In this work we presented one such simulation, imple¬ 
menting repetitive QEC with the [[5,1,3]] code on a six- 
qubit network with always-on Ising couplings and clas¬ 
sical correlated phase noise as a source of decoherence. 
The one- and two-qubit gates were implemented via care¬ 
fully designed sequences of shaped pulses. Realistically 
simulating such gates and associated errors requires inte¬ 
grating the corresponding multi-qubit unitary dynamics. 
Such simulations, like current experiments, are limited to 
very small system sizes. As a result, one can use only sim¬ 
plest weak codes, with very few ancillary qubits, which 
puts additional constraints on the accuracy of the imple¬ 
mented gates. 

As in any case where gates are designed perturba- 
tively, up to some fixed order in the perturbation (in¬ 
teraction) Hamiltonian, the systematic errors associated 
with our gates are correlated multi-qubit errors, which 
worsen with the increased connectivity of the qubit net¬ 
work. On the other hand, the [[5,1,3]] code we used is 
able to correct only single-qubit errors. To make QEC 
possible, we had to tune the couplings down and make 
the two-qubit gates longer, increasing the intrinsic fideli¬ 
ties of our gates to six nines or more. As a result, just a 
few rounds of repetitive QEC required tens of thousands 
of pulses, with the fidelity noticeably suffering, e.g., from 
relatively modest integration errors (not shown). In this 
weakly coupled regime, our simulations show that single¬ 
qubit errors due to phase noise do not propagate exces¬ 
sively. 

Overall, we demonstrated repetitive quantum error 
correction in a fully quantum-mechanical simulation, 
with the error correction responsible for the average in¬ 
fidelity reduction by an order of magnitude or more. We 


have also presented a combined protocol integrating dy¬ 
namical decoupling and quantum error correction. Dy¬ 
namical decoupling is particularly effective against low- 
frequency noise which in our simulations had an asymp¬ 
totic dephasing time T 2 as short as few nominal pulse 
lengths Tp. We also saw that our combined DD/QEC 
protocol remains effective in the presence of a weak high- 
frequency phase noise. 

While dephasing-only model appear to be too simplis¬ 
tic, we notice that as a result of controlled dynamics, 
some of the dephasing propagates to the longitudinal 
channel |46j. In particular, our original simulations which 
involved similar gates with the three-qubit code protect¬ 
ing against single-qubit phase errors, produced a much 
smaller fidelity improvement due to QEC [73]. 

Our model excludes many physical effects which may 
be relevant for engineering a quantum computer based on 
a specific qubits implementation, such as couplings be¬ 
tween nominally disconnected qubits, multi-level struc¬ 
ture of the solid-state qubits and corresponding leakage 
errors, violations of the rotating wave approximation, re¬ 
alistic decoherence which may produce additional cor¬ 
relations between the qubits, etc. Even when the cor¬ 
responding effects are small, they can result in errors 
correlated in time or between qubits, and thus strongly 
affect the overall coherent multi-qubit dynamics. De¬ 
signing coherence protection schemes with improved sta¬ 
bility against such effects is also possible, if one knows 
which decoherence mechanisms are dominant. Each ad¬ 
ditional improvement would require more finely tuned 
pulse shapes, longer gates, or a longer code, increasing 
the requirements on the dynamical range of the qubit 
system used in the experiment. Thus, in our opinion, 
careful studies of realistic models which incorporate such 
effects are absolutely necessary in order to construct a 
scalable quantum computer. 
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